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Abstract 

We complete the existing literature on the structure and stability of polytropic gas 
spheres reported in the classical monograph of Chandrasekhar (1942). For isolated poly- 
tropes with index 1 < n < 5, we provide a new, alternative, proof that the onset of 
instability occurs for n = 3 and we express the perturbation profiles of density and veloc- 
ity at the point of marginal stability in terms of the Milne variables. Then, we consider 
the case of polytropes confined within a box of radius R (an extension of the Antonov 
problem for isothermal gas spheres). For n > 3, the mass-density relation presents some 
damped oscillations and there exists a limiting mass above which no hydrostatic equilib- 
rium is possible. Like for isothermal gas spheres, the onset of instability occurs precisely 
at the point of maximum mass in the series of equilibrium. Analytical results are obtained 
for the particular index n = 5. We also discuss the relation of our study with generalized 
thermodynamics (Tsallis entropy) recently investigated by Taruya & Sakagami (cond- 
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1 Introduction 

In earlier papers (Chavanis 2001a,b) we investigated the gravitational instability of finite 
isothermal spheres in Newtonian gravity and general relativity. On a theoretical point of view, 
these systems exhibit an interesting behavior with the occurence of phase transitions associated 
with gravitational collapse. A rich stability analysis follows and can be conducted analytically 
or by using graphical constructions. On an astrophysical point of view, these studies can be rel- 
evant for various systems including elliptical galaxies, globular clusters, the interstellar medium 
or the core of neutron stars. In this paper, we propose to extend our study to the case of poly- 
tropic gas spheres. Polytropes with index 1 < < 5 are self-confined, so it is not necessary to 
introduce an artifical "box" to limitate their spatial extent. Using the methods developed for 
isothermal configurations, we show that the transition from stability to instability corresponds 
to an index n = 3. This result is well-known but we provide a new derivation based on the 
exact resolution of the pulsation equation for polytropes. The perturbation profiles at the point 
of marginal stability are expressed in terms of the Milne variables. The density profile has only 
one node and the velocity perturbation is proportional to the radial distance. Then, we consider 
the case of polytropes with arbitrary index n > 1 confined within a box of radius R. For n = oo, 
we recover the classical Antonov (1962) problem for isothermal gas spheres. However, for n > 3 
we already obtain results strikingly similar to those obtained for isothermal configurations. In 
particular, confined polytropes in hydrostatic equilibrium can exist only below a limiting mass 
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(for a given box radius R) and the series of equilibrium becomes unstable precisely at the point 
of maximum mass. Subsequent oscillations in the mass-density profile (for n > 5) are associated 
with secondary modes of instability. The locii of these modes of instability follow a geometric 
progression with a ratio depending on the index of the polytrope. For n = oo, we recover the 
ratio 10.74... of isothermal gas spheres. While this paper was in preparation, we came accross 
the preprint of Taruya & Sakagami (2001) on a related subject. These authors investigate the 
stability of polytropes in the framework of extended thermodynamics, using Tsallis entropy. 
Therefore, in section ^, we analyze the connexion of the present study with their approach and 
we discuss the relevance of Tsallis entropy for describing astrophysical systems. 

2 Properties of polytropic gas spheres 
2.1 The Lane-Emden equation 

Polytropic stars are characterized by an equation of state of the form 

(1) p = Kp\ 

where K and 7 are constants. The index n of the polytrope is defined by the relation 

(2) 7 = 1 + -- 

n 

The equation of state (^ corresponds to an adiabatic equilibrium in regions where convection 
keeps the star stirred up and produces a uniform entropy distribution (s = Cst.). In that case, 
7 is the ratio of specific heats Cp/cy at constant pressure and volume. For a monoatomic gas, 
7 = 5/3. The equation of state (|l|) also describes a polytropic equilibrium characterized by 
a uniform specific heat c = 6Q/dT. In this more general situation 7 = (cp — c)/(cy — c). 
Convective equilibrium is recovered for c = and isothermal equilibrium is obtained in the 
limit of infinite specific heat c — > +00. A power law relation between pressure and density 
is also the limiting form of the equation of state describing a gas of cold degenerate fermions 
(Chandrasekhar 1942). In that case, the constant K can be expressed in terms of fundamental 
constants. In the classical limit 7 = 5/3, n = 3/2 and K = ^(3/7r)^/^/i^/m (where h is the 
Planck constant) and in the relativistic limit 7 = 4/3, n = 3 and K = |(3/7r)^/^/ic (where 
c is the speed of light). Historically, the index 5/3 appears in the classical theory of white 
dwarf stars initiated by Fowler (1926) and the index 4/3 is related to the limiting mass of 
Chandrasekhar (1931). 

The condition of hydrostatic equilibrium of a spherically distribution of matter reads 

(3) ^ = -p— 

dr dr ' 

where $ is the gravitational potential. Using the Gauss theorem 

d^ GM{r) 



(4) 



dr r 
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where M(r) = pAnr'^dr is the mass contained within the sphere of radius r, we can derive 
the fundamental equation of equilibrium (Chandrasekhar 1942) 

(5) Id/ dp\ _ ^^^^ 

r"^ dr\ p dr J 
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Eqs. (P® fully determine the structure of polytropic gas spheres. Letting 



(6) 



P = PoO'^ 



l-l/n 



1/2 



_ K{n + 1) J 

where po is the central density, we can reduce the condition of hydrostatic equilibrium to the 
Lane-Emden equation (Chandrasekhar 1942) 



(7) 



with boundary conditions 

(8) ^(0) = 1, ^'(0) = 0. 

For n > 3, the Lane-Emden equation admits an analytical solution which is singular at the 
origin: 



(9) 



2(n-3) 
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Regular solutions of the Lane-Emden equation must in general be computed numerically. For 
^ — s> 0, we can expand the function 9 in Taylor series. The first terms in this expansion are 
given by 



(10) 



6^ 120^ 



The behavior of 9{^) at large distances deserves a more specific discussion. For 1 < n < 5, the 
density falls off to zero at a finite radius i?, identified as the radius of the star. We denote by 
^1 the value of the normalized distance at which ^ = 0. For ^ ^i, we have 
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For n > 5, the polytropes extend to infinity, like the isothermal configurations recovered in the 
limit n = oo. For ^ +oo, 
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The curve (^2|) intersects the singular solution infinitely often at points that asymptotically 
increase geometrically in the ratio 1 : exp{2(n — 1)ti/ \J1v? — 22n — 1}. Since 9"- ~ at 
large distances, these configurations have an "infinite mass", which is clearly unphysical. In 
the following, we shall restrict these configurations to a "box" of radius R, like in the classical 
Antonov problem. Therefore, Eq. (|^ must be solved for ^ < a with 

L A (n + 1) 

Note that, for a fixed box radius i?, a is a measure of the central density po- The case n = 5 
is special. For this index, the Lane-Emden equation can be solved analytically and yields the 
result: 

The total mass of this configuration is finite but its potential energy diverges. Therefore, this 
polytrope must also be confined within a box. On Fig. |l|, we have ploted different density 
profiles corresponding to the polytropic indices n = 3, 5, 6. 
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Figure 1: Density profiles of polytropes with index n = 3,5, 6. The dash fine corresponds to 
the singular solution @. 

2.2 The Milne variables 

Like in the analysis of isothermal gas spheres, it will be convenient in the following to introduce 
the Milne variables {u,v) defined by (Chandrasekhar 1942): 

.... ^0' 

(15) « = -^, v = --. 

Taking the logarithmic derivative of u and v with respect to ^ and using Eq. (|^), we get 

(16) = h3-nv -u). 

Due to the homology invariance of the polytropic configuations (see Chandrasekhar 1942), the 
Milne variables u and v satisfy a first order differential equation 

, , u dv u + V — 1 

(18) — — = . 

V du u + nv — 3 

The solution curve in the {u, v) plane (see Figs. Q-^) is parametrized by ^. It starts from the 
point {u,v) = (3,0) with a slope {dv/du)o = —-^ as ^ 0. For 1 < n < 5, the curve is 
monotonous and tends to {u,v) = (0, +oo) as ^ — > ^i. More precisely, using Eq. (|llD, we have 

71+ 1 

(19) uv'^r^u:-', u;n = -^r(^'i 

For n > 5, the solution curve spirals indefinitely around the fixed point {Us,Vs) = {^E^, ;7ri)) 
corresponding to the singular solution (|]), as ^ tends to infinity. All polytropic spheres must 
necessarily lie on this curve. For bounded polytropic spheres, ^ must be terminated at the box 
radius a. For n = 5, we have the following analytical relation between the Milne variables: 

(20) u + 3v = 3. 
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Figure 2: The {u,v) plane for polytropes with index 1 < n < 5. The construction is made 
exphcitly for n — A. 




Figure 3: The {u,v) plane for polytropes with index n = 5. 
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Figure 4: The (m, v) plane for polytropes with index n > 5. The construction is made exphcitly 
for n = 6. 

2.3 The maximum mass and minimum temperature of confined poly- 
tropes 

For polytropes confined within a box of radius i?, there exists a well-defined relation between 
the mass M of the configuration and the central density po (through the parameter a). Starting 
from the relation 



(21) 



R 

M = I pinr^dr = 47rpo 

'0 



l/n-l 




AttG 



3/2 



and using the Lane-Emden equation (|^), we get 



(22) 



M = -inpo 



AtiG 



Expressing the central density in terms of a, using Eq. (p!3[), we obtain after some rearrange- 
ments 



K{l+n) 



(23) M = -Att 

Introducing the parameter 

(24) 

the foregoing relation can be rewritten 

(25) 7] = -a^^'(a). 
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For n < 5, the normalized box radius a in necessarily restricted by the inequality a < ^i. For 
the limiting value a = ^i, corresponding to an isolated polytrope satisfying p{R) = 0, we have 

(26) 77(6) = cu„. 

The quantity Un, defined by Eq. ([19|), is tabulated in Chandrasekhar (1942). The definition 
( p^ ) of f] and the relation ( p5| ) between 77 and a are consistent with the formulae derived in 
the case of an isothermal gas (see Chavanis 2001a). This connexion is particularly relevant if 
we interpret the constant as a polytropic temperature 6^ (see Chandrasekhar 1942, p. 86). 
For n — ^ +00, we find that rj ~ rj^/n where rj^ = 1^^^^ is the normalized temperature of an 
isothermal sphere (/? = 1/kT). 
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Figure 5: Mass-density profiles for polytropic configurations with index n = 2, 3, 4, 5, 6. A mass 
peak appears for the first time for the critical index n = 3. 

The function r]{a) is represented on Fig. ^ for different values of the polytropic index n. 
Instead of a, we could have used the density contrast 

(27) n^^^^Ha)-". 

Using the Lane-Emden equation (^ and the definition of the Milne variables (0), it is straight- 
forward to check that the condition of extremum drj / da = is equivalent to 

n — 3 

(28) Mo = = Us, 

n — 1 

where, by definition, uq = u{a) and Ug refers to the singular solution (H). For n — 00, we 
recover the condition uq = 1 previously derived for isothermal configurations (Chavanis 2001a). 
The values of a for which rj is extremum are given by the intersection(s) between the solution 
curve in the (u, v) plane and the line u = Ug (see Figs. @-§). For 1 < n < 3, Us < and there is 
no intersection. The mass-density relation is therefore monotonous. For 3 < n < 5, the curve 
r){a) presents a single maximum at ai. For n = 3, this maximum is reached at the extremity 
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of the curve (ai = ^i). For n = 5, ^ +00 and the function 77(a) is exphcitely given by 



(29) 



V 



3(l + ia2)3/2' 



The maximum of rj is located at ai = yl5. Finally, for n > 5, the mass-density relation 
presents an infinite number of damped oscillations since the line defined by Eq. (|28| ) passes 
through the center of the spiral in the {u, v) plane. If we denote by ak the loci of the extrema 
of rj{a), these values asymptotically follow the geometric progression 



(30) 



2(n-l)7r 
) \/7>i2-22n- 



{k — > 00, integer). 



obtained by substituting the asymptotic expansion (|T^ in Eq. (P5|). For n ^ 00, we recover 
the ratio e^'^l^ = 10.74... corresponding to classical isothermal configurations (Semelin et al. 
1999, Chavanis 2001a). 

From the above results, it is clear that restricted polytropic spheres with index > 3 can 
exist only for 



(31) 



n+l 



This implies in particular the existence of a limiting mass (for a given confining radius R) such 
that 



(32) 



M < Mrr 



47T7]„ 



K{l + n) 
4nG 



71 — 3 

R—. 



Alternatively, for a given mass M and radius R, this inequality implies the existence of a 
minimum value of the polytropic temperature 0^ = Km/k. Indeed, 0^ is restricted by the 
inequality 



(33) 



07 > 



ATcGm 



k(l + n) xAttt]. 



M 



1 



R' 



In the limit n = 00, we recover the critical temperature kTc 



GMm 
2.52R 



below which an isother- 



mal sphere is expected to collapse (Lynden-Bell & Wood 1968). This minimum temperature, 
corresponding to drj/da = 0, appears for the first time for a critical index n = 3. This obser- 
vation will take a deeper physical significance in the stability analysis of the following section. 
For 77. < 3, the mass-density profile is monotonous and t] < ujn- The maximum mass (resp. 
minimum temperature) of confined polytropes is always smaller (resp. larger) than the corre- 
sponding one for isolated polytropes. However, this bound does not correspond to a condition 
of extremum drj/da = but rather to the impossibility of constructing polytropes with a > ^1 
(since 6 can become negative). 



3 Dynamical stability of polytropic gas spheres 
3.1 The equation of radial pulsations 

We shall now investigate the dynamical stability of polytropic gas spheres against radial per- 
turbations. We use a method and presentation similar to that adopted for isothermal spheres 
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in Chavanis (2001a). The equations describing the motions of a gaseous star are the equation 
of continuity, the Euler equation and the Poisson equation 

(34) ^ + V(pv) = 0, 

dv 1 

(35) — + (vV)v = — Vp - V$, 

ot p 

(36) A* = AnGp. 

It is possible to show that viscosity does not change the onset of instabihty. Therefore, for 
simphcity, we have directly used the Euler equation instead of the Navier-Stokes equation. We 
assume furthermore that the pressure and the density are connected by the polytropic equation 
of state d^). We write small perturbations around equilibrium in the form 

(37) v = (5v(r,t), p = p + V(r,t), 



(38) p = p + 5p{r,t), $ = $ + (5<l>(r,t), 

where the bar refers to the stationary solution (in the following we shall drop the bar). The 
linearized equations for the perturbations are 

(39) p-^ = -K-iV{p^-^5p) - 5pV$ - pV(5$, 

(40) ^ + V(p5v) = 0, 

(41) A^^ = A-nG5p. 

Restricting ourselves to radial perturbations and writing the time dependance of the perturba- 
tion in the form 5v ~ e^*, 5p ~ e''**,..., the equations of the problem become 

(42) Xp5v = -K^^y-^5p) - 6p^ - p^, 

(43) X5p+\^{pr^5v) = 0, 

44 U2 =4^G'5p. 

dr \ dr J 

We introduce the function q{r) by the relation 

47rr^ dr 
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Physically, g(r) represents the mass perturbation SM{r) = SpAirr'^dr within the sphere of 
radius r. Thus, by definition, g(0) = 0. Then, it is readily seen that the Poisson equation 
is equivalent to 

(46) '-^ = % 

ar 

which is just the perturbed Gauss theorem (^). On the other hand, the continuity equation 
(143|) leads to the relation 



(47) 6v = --^g. 

Substituting these results back into Eq. (|42|), we obtain 

, d f ^ . I dq\ 1 dqd^ Go 



47rr^ dr \ ATcr"^ dr J Anr'^ dr dr 

Using the condition of hydrostatic equilibrium (|^), the foregoing equation can be rewritten 

±fp^dq\ Gq _ 



(49) ^7:r + ^ = i — 2i- 

dr \47rr^ dr J r^ Airpr"^ 
which is the required equation of radial pulsations for a polytrope. 

3.2 The condition of marginal stability 

Considering the case of marginal stability A = and introducing the dimensionless variables 
defined in section H, we can reduce the equation of radial pulsations to the form 

d fe^-''dF\ nF ^ 

with -F(O) = 0. Denoting by C the differential operator appearing in Eq. (|50|) and using the 
Lane-Emden equation (|^, we easily establish that 

rl rl \ rl 



(52) cieon = ^ i^-^^ieon) + = (3 + n)e' + ine" + nie^ = (3 - n)e'. 

Therefore, the general solution of Eq. (|50|) is 

(53) F(o = c, (^ee- + "^/e^ , 

where Ci is an arbitrary constant. The connexion with isothermal configurations (see Chavanis 
2001a) is particularly obvious if we make the correspondance 9' ip' and ^" e~'^. 
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3.3 Boundary conditions 

The equation of pulsations ( ^9|) must be supplemented by appropriate boundary conditions. 
For self-confined polytropes (1 < n < 5), we require that the Lagrangian derivative of the 
pressure vanishes at the surface of the configuration, i.e. 

(54) ^6p + 6v^ = 0, atr = R. 

at dr 



(55) ^(5-^^1-0. atr^fl. 



Using Eqs. (^) (|47|) , this condition can be rewritten 

dq q dp 
dr p dr 

or, in dimensionless form, 

dF d9 

(56) e- - n-F = 0. aU = 6- 

The index of the marginally stable polytrope is obtained by substituting the solution (|53|) in 
the boundary condition (|56D . Since 6{C,i) = 0, this yields n = 3. Therefore, the transition from 
stability to instability occurs for a polytropic index n = 3 or, equivalently, for an adiabatic index 
7 = 4/3. Of course, this result is well-known and could have been obtained from the general 
theorems of stellar pulsation (see, e.g.. Cox 1980). However, we provide here an alternative 
method based on the explicit resolution of the equation of pulsations for polytropes. 

According to Eq. (^5]), the perturbation profile at the point of marginal stability is given 

by 

(57) ^ = 



with the expression (|53D for F{^). Simplifying the derivative with the aid of the Lane-Emden 
equation (|^), we obtain 



where v is the Milne variable (|T5D. For n = 3, it reduces to 

(59) Sp=^cipoe'il-v). 

47r 

This perturbation profile is ploted on Fig. On the other hand, the velocity profile is given 
by Eq. (|47|). Introducing the velocity of sound c^, we find that 



^ ^ CS 4:71 - 1 

For n = 3, we get 

(61) — = -^ci^, 

cs 47r 

and we observe that the velocity is proportional to the radial distance. 
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Figure 6: Density perturbation profile for the critical polytrope of index n = 3. 



For polytropes confined within a box, the boundary condition consistent with the conser- 
vation of mass is F{a) = 0. With the expression ( [53| ) for -F(^), this yields 

3 



(62) 



n 



Uo 



n 



which is precisely Eq. (pSl). Therefore, the point of marginal stability coincides with the point 
of maximum mass (or minimum temperature) in the series of equilibrium. For n < 3, restricted 
polytropes are always stable since the function rj{a) is monotonous. For 3 < n < 5, the 
configurations with a > ai are unstable but there are no secondary instabilities. For n > 5, a 
first mode of stability is lost for a = ai. Subsequent oscillations in the mass-density profile are 
associated with secondary instabilities. 

The profiles of density and velocity that trigger these instabilities can be written 



(63) 



5p 
P 



v), 



(64) 



6v 
Cs 



y 9' , 

47r Q'^^ 



u 



where we recall that ^ < ak for the mode of order k. We can determine the qualitative behavior 
of these profiles graphically by considering the intersections between the solution curve in the 
{u,v) plane and the lines v = Vg and u = Ug (see Figs. 0-^. These intersections determine the 
nodes of the perturbation profile. For n < 5, there is only one intersection since the {u, v) curve 
is monotonous. Therefore, the density profile 6p/p has only one node at .^(i) = 2.16... (see Fig. 
For n = 5, we obtain the analytical solutions 



(65) 



~p ~ 8^^^ 3 + ^2' 



(66) 
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The perturbation 5p becomes zero at ^(i) = y^- For n > 5, there are several intersections 
with the spiral since the lines v = Vs and u = Ug pass through the center of the spiral. The 
description of the perturbation profiles is similar to the one given for isothermal spheres. In 
particular, the fundamental mode of instability has only one node and high order modes of 
instability present numerous oscillations whose nodes follow a geometric progression. We refer 
to our previous paper (Chavanis 2001a) for a more precise description of these results. 

4 Generalized thermodynamics and Tsallis entropy 

Recently, Taruya & Sakagami (2001) have investigated the stability of polytropic spheres within 
the framework of generalized thermodynamics. It has been argued by Tsallis and co-workers 
that ordinary statistical mechanics and thermodynamics does not describe correctly systems 
with long-range interactions, which are by essence non extensive. A family of functionals of the 
form 

(67) S, = j ir - f)dhd\ 

known as Tsallis entropies, has been introduced to extend the classical Boltzmann-Gibbs sta- 
tistical mechanics to these systems. These functionals are labeled by a parameter q. The 
Boltzmann entropy is recovered in the limit g — 1. This new formalism has been applied in 
various domains of physics, astrophysics, fluid mechanics, biology, economy etc... and it extends 
the results obtained with ordinary statistical mechanics. In the case of self-gravitating systems, 
it was shown by Plastino & Plastino (1993) that the extremalization of the Tsallis entropy at 
fixed mass and energy yields a polytropic equation of state of the form (|l|). The parameter q 
is related to the index n of the polytrope by the relation 

(68) n = ^ + |, (g>l). 

q — 1 2 

Then, Taruya & Sakagami (2001) investigated the stability problem in the microcanonical en- 
semble by extending the analysis of Padmanabhan (1989) for the classical Antonov instability. 
Polytropic configurations are said to be stable if they correspond to (generalized) entropy max- 
ima at fixed mass M and energy E. Taruya & Sakagami showed that, for n > 5, an equilibrium 
exists only above a critical energy depending on the index of the configuration. In the limit 
n — > +00, the Antonov result = —0.335GM^/R obtained with the Boltzmann entropy is 
recovered. Furthermore, they showed that the onset of instability coincides with the point of 
minimum energy in the series of equilibrium in agreement with standard turning point analysis. 
For n < 5, there is no critical value of energy and the polytropic configurations are stable. These 
results differ from those obtained in the present paper where it is found that the index marking 
the transition from stability to instability is n = 3 (in agreement with standard theorems of 
stellar pulsations). The origin of this discrepency is certainly related to the inequivalence of 
the statistical ensembles in extended thermodynamics, like in ordinary thermodynamics, for 
self-gravitating systems (see, e.g., Padmanabhan 1990). Taruya & Sakagami work in the mi- 
crocanonical ensemble and maximize Tsallis entropy at fixed mass and energy. On the contrary, 
in our study, we keep the polytropic temperature = Km/k constant, which corresponds 
to the canonical description. It would be of interest to study the maximization of Tsallis free 
energy at fixed mass and temperature to see if it provides the same conditions of stability as 
our dynamical approach based on the Navier-Stokes equations. For the Boltzmann entropy, a 
clear connexion was found between dynamical stability and thermodynamical stability in the 
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canonical ensemble (Semelin et al. 2001, Chavanis 2001a) and it is desirable to check whether 
this connexion is preserved by Tsallis generalized thermodynamics. This requires in particular 
to define properly the notion of temperature and free energy in the generalized sense. 

On Figs. [7|-PH|, we have represented the equilibrium phase diagram (for different values of 
n) obtained in the framework of extended thermodynamics, interpreting the parameter 77 as a 
generalized inverse temperature. The temperature-energy curve is defined in a parametric form 
by the equations 



(69) 



n + l 

-a"-^6 [a] 



(70) 



A 



ER 
GAP 



n 



3/ I \ n-2u{a) 
2 \ v{a) J ^ n + l v{a) 



The expression ([70|) for the energy has been derived by Taruya & Sakagami (2001). For n < 3, 
there is no turning point in the diagram, so the polytropes are always stable. For 3 < n < 5, the 
inverse temperature t] presents a maximum but not the energy. Therefore, the polytropes are 
always stable in the microcanonical ensemble but they are unstable in the canonical ensemble 
after the turning point. It can be noted that the unstable region has a negative specific heat 
C ~ ^ < (in the generalized sense). For n > 5, the temperature and the energy both present 
an infinite number of extrema. For n ~ 6, there are crossing points at which two solutions have 
the same value of temperature and energy (but a different density contrast). The polytropes 
are unstable in the canonical ensemble after the first turning point of temperature and they are 
unstable in the microcanonical ensemble after the first turning point of energy. These results 
are strikingly similar to those obtained with the classical Boltzmann entropy {n = 00). 



n=2 
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Figure 7: Phase diagram for n = 2. 



It is interesting that the results obtained for isothermal gas spheres described by ordinary 
thermodynamics (Boltzmann entropy) can be extended to polytropic configurations by intro- 
ducing a generalized functional (Tsallis entropy). This makes this new formalism of interest. 
However, we would like to make two comments, (i) First of all, Tsallis entropy does not describe 
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Figure 10: Phase diagram for n = 10. 



polytropic gaseous stars although it yields an equation of state of the polytropic form. Indeed, 
Tsallis entropy predicts in addition a distribution function of the form 



whereas in polytropic stars, the distribution function corresponds to a local thermodynamical 
equilibrium (L.T.E) characterized by a MaxweUian distribution of velocities 



with a local temperature T(r) (recall that, for convective equilibrium, the entropy is uniform 
not the temperature). In view of these remarks, the connexion between (generalized) ther- 
modynamical stability (in the canonical ensemble) and dynamical stability with respect to 
Navier-Stokes equations is not obvious. It is however expected to be true since the condition of 
dynamical stability is consistent with the turning point argument of thermodynamical stability, 
(ii) Tsallis entropy could describe a particular class of stellar systems, apparently not observed 
in nature, known as stellar polytropes (see Binney & Tremaine 1987). These systems have the 
same structure as polytropic gas spheres in physical space but not in phase space (the equiv- 
alence only occurs for a polytropic index n = oo, i.e. for an isothermal configuration). This 
distinction is important because isolated stellar polytropes are always stable with respect to 
the Vlasov (or coUisionless Boltzmann) equation (see Binney & Tremaine 1987) while isolated 
polytropic stars are unstable with respect to the Navier-Stokes equations for n > 3. General- 
ized thermodynamics could be developed (at least formally) for this particular class of stellar 
systems. 

However, generalized thermodynamics is often presented as a way of solving the "problems" 
associated with the use of the Boltzmann entropy in systems with long-range interactions and 
we would like to criticize this interpretation. For self-gravitating systems, the "problems" 
usually reported are the absence of global entropy maximum (resulting in the "infinite mass 



(71) 




(72) 
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problem" and the "gravothermal catastrophe"), the negative specific heats, the inequivalence of 
statistical ensembles and the absence of thermodynamical limit (a discussion of these concepts 
can be found in Chavanis et al. 2001). It should be noted that the Tsallis entropy displays 
the same phenomena, at least for values of g < 9/7 (i.e. n > 5) considered by Taruya & 
Sakagami (2001). On the other hand, these relatively unusual phenomena should not cause 
surprise if one recognizes that self-gravitating systems are relatively exceptional among A^- 
body systems. For example the "infinite mass problem" (see, e.g., Binney & Tremaine 1987) 
is just a mathematical curiosity. There is no justification in maximizing the entropy in an 
infinite domain even if an entropy maximum exists. Due to kinetic effects, the relaxation 
is necessarily incomplete and the ergodic hypothesis which sustains the statistical analysis is 
restricted to a finite region of space. This incomplete relaxation is qualitatively discussed by 
Lynden-Bell (1967) in his statistical description of the "violent relaxation" of collisionless stellar 
systems (e.g., elliptical galaxies). A similar limitation occurs in the context of two-dimensional 
turbulence to understand the confinement of vortices that form after a rapid merging (see, e.g., 
Chavanis & Sommeria 1998 and in particular Brands et al. 1999 for a discussion of Tsallis 
entropy). This incomplete relaxation can be described by kinetic equations of the form given 
by Chavanis et al. (1996) with a space dependant diffusion coefficient. A convincing numerical 
evidence of this kinetic confinement was given by Robert & Rosier (1997) in two-dimensional 
turbulence. In addition, in astrophysics, specific processes (e.g., the evaporation of stars, tidal 
effects,...) must be taken into account and can restrict the domain of applicability of statistical 
mechanics. These limitations are usually accounted for by introducing truncated models like 
the Michie-King model for globular clusters or the models proposed by Stiavelli & Bertin 
(1987) or Hjorth & Madsen (1993) for elliptical galaxies. These models lead to composite 
configurations with an isothermal core and a polytropic envelope. With these modifications 
(justified by precise physical arguments), ordinary statistical mechanics provides in general a 
good explanation of astrophysical phenomena and is consistent with observations. Note that 
the above-mentioned truncated models differ from pure polytropes, so that they cannot be 
justified by generalized thermodynamics Q. On the other hand, the "gravothermal catastrophe" 
(Lynden-Bell & Wood 1968) should not throw doubts on the validity of statistical mechanical 
arguments applied to self-gravitating systems. On the contrary, it is gratifying that standard 
thermodynamics is consistent with the natural tendency of self-gravitating systems to undergo 
a gravitational collapse (Jeans instability). The gravothermal catastrophe has been confirmed 
by sophisticated numerical simulations (see, e.g., Larson 1970, Cohn 1980) based on standard 
kinetic theory (Landau-Fokker-Planck equations) and it is expected to be at work in globular 
clusters. The existence of negative specific heats in the microcanonical ensemble (and the 
resulting inequivalence of statistical ensembles) is a direct consequence of the Virial theorem for 
self-gravitating systems and is not in contradiction with ordinary thermodynamics (Lynden-Bell 

^Recently, Fa & Pedron (2001) have proposed to describe elliptical galaxies by an extended King model 
using Tsallis entropy. They obtained a good fit for NGC 3379. However, we beleive that this agreement is 
fortuitous. An even better fit has been obtained by Hjorth & Madsen (1993) with a different model justified by 
more convincing physical arguments. Their model has an isothermal (instead of a polytropic) inner region in 
agreement with Lynden-Bell's statistical mechanics of violent relaxation. The polytropic envelope is justified by 
incomplete relaxation (or equivalently by a lack of ergodicity) when we depart from the central mixing region. 
Their model differs from the Michie-King model because the system is collisionless and not subject to tidal 
forces (King's models, which can be derived rigorously from the Fokker-Planck equation, only apply to tidally 
truncated coUisional stellar systems such as globular clusters). It is seen on Fig. 3 of Fa & Pedron (2001) that 
an isothermal distribution in the inner region provides a sensibly better fit to observations than a polytropic 
one (it can reproduce the oscillations). In addition, the model of Hjorth & Madsen (1993) has only one fitting 
parameter $o (the central potential) instead of two $0 and g in Fa & Pedron's analysis. These remarks tend to 
show that generalized thermodynamics is not necessary to understand the structure of self-gravitating systems. 
A similar conclusion has been reached by Brands et al. (1999) in two-dimensional turbulence. 
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& Lynden-Bell 1977, Padmanabhan 1990). Finally, self-gravitating systems have a well-defined 
(albeit unusual) thermodynamical limit in which the number of particles N and the volume 
go to infinity keeping N/R fixed (de Vega & Sanchez 2001). 

Therefore, we are tempted to beleive that ordinary statistical mechanics is still relevant for 
describing self-gravitating systems and two-dimensional vortices. An isothermal distribution 
corresponds to the most probable distribution reached by a system after a complex evolution 
during which microscopic information is lost. Accordingly, it maximizes the Boltzmann en- 
tropy, which is a measure of the number of microstates associated with a given macrostate (see, 
e.g., Ogorodnikov 1965, Lynden-Bell 1967). This definition of the entropy does not depend 
whether the particles are in interaction or not. Nonlocality intervenes only when the entropy 
is maximized for a given value of energy. It is implicitly assumed, however, that all accessible 
microstates are equiprobable, which is the fundamental postulate of statistical mechanics. Gen- 
eralized thermodynamics is interesting for extending standard results of thermodynamics to a 
larger class of functionals which possess nice mathematical properties, but it should not be pre- 
sented (in our point of view) as an improvement or a replacement of classical thermodynamics. 
Wc agree with Tsallis and co-workers that the kinetic theory of self-gravitating systems (see 
Kandrup 1981) and two-dimensional vortices (see Chavanis 2001c) encounters some problems 
due to the occurence of memory effects, spatial delocalization and logarithmic divergences in 
the diffusion coefficient. For these reasons, and also for the problems of incomplete relaxation 
(lack of ergodicity) mentioned previously, the system may not necessarily reach the maximum 
entropy state described by the Boltzmann distribution, even in the mixing region. The de- 
viation from the Boltzmann distribution is never too severe (see, e.g., the isothermal core of 
globular clusters and of elliptical galaxies) so that ordinary statistical mechanics provides a 
fairly good prediction from any initial condition. It may happen that the effective equilibrium 
distribution can be fitted by the g-distribution proposed by Tsallis and co-workers better than 
by the Boltzmann distribution. This is the case for some examples of statistical equilibrium in 
two-dimensional turbulence (see the discussion of Brands et al. 1999). However, this agreement 
may be fortuitous rather than dictated by a general physical principle since there is a parameter 
q in the theory which, in practice, must be adjusted in each case. In fact, it is advocated by 
Tsallis (private communication) that this parameter is not free but uniquely determined by 
the microscopic dynamics of the system. It can be relatively easy to determine if the system 
is simple enough but it can also be very hard to determine in the case of complex systems like 
those involved in fluid mechanics and stellar dynamics. In that case, it must be regarded as 
a fitting parameter. Unfortunately, the power of prediction of Tsallis theory is limited as long 
as a precise prescription for determining q is not given. Tsallis entropy could be considered, 
however, as a heuristic attempt to take into account non ergodic effects in complex systems. In 
view of these remarks, it would be of considerable importance to see if the non Markovian and 
non local generalized kinetic equations proposed by Kandrup (1981) in stellar dynamics and by 
Chavanis (2001c) in vortex dynamics allow for a stationary distribution of the form conjectured 
by Tsallis, instead of the Boltzmann distribution. From our point of view, we beleive that they 
do not select any "universal" distribution except in the approximation where (i) memory ef- 
fects arc ignored (ii) a local approximation is made (in the stellar context) (iii) ergodicity is 
assumed, in which case the classical Boltzmann distribution is obtained. It is hard to beleive 
that extremely complicated effects of non ergodicity, spatial delocalization and memory can be 
encapsulated in a simple functional with a prescribed value of q. However, if it can be proved 
that these generalized kinetic equations rigorously converge towards a TsaUis distribution with 
g 7^ 1, this would of course be a strong argument in favour of generalized thermodynamics. 
Unfortunately, the study of these generalized kinetic equations seems of considerable difficulty 
and will demand a long standing effort. 



18 



5 Conclusion 



This study has revealed that the striking results obtained by Antonov (1962) and Lynden-Bell 
& Wood (1968) in their investigations on the thermodynamics of self-gravitating systems are 
not restricted to isothermal configurations but also occur for polytropic configurations. It is 
possible that this formal analogy be the mark of a generalized thermodynamics (Tsallis entropy) 
but this point needs to be discussed in greater detail. In any case, the study presented in this 
paper and in our companion papers (Chavanis 2001a,b) provide a unified description of the 
structure and stability of isothermal and polytropic gas spheres. This is a useful complement 
to the existing literature on the subject (Chandrasekhar 1942). Quite remarkably, the stability 
analysis can be performed analytically or by using graphical constructions. Our study is limited 
in its applications by the requirement of a confining box in the case of isothermal configurations 
and polytropes of index n > 5. However, this box model (first introduced by Antonov) should 
be sufficient to exhibit the essential features of these systems and it is also very convenient for a 
theoretical analysis. In future works (in preparation), we shall extend our analytical results to 
the case of composite configurations made of an isothermal core and a polytropic envelope with 
index n < 5. This is an important generalization for a better description of self-gravitating 
systems from the viewpoint of statistical mechanics and thermodynamics. 
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